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ABSTRACT 


The NPS computer simulation model was modified to study 
sputtering by molecular ions. The simulations were 
performed on a Cu(11ll) surface for normally incident 05 
NOS O AS Sana 5.0 kev/ion. The molecule's 
angular orientation was both set at specific values and 
randomly determined. The normalized energy distributions 
of sputtered atoms for two energy ranges, 0-20 and 0-3 
ae spubtesrea atom mand the sputtering yield ratios were 
investigated. The 0-20 ev/sputtered atom energy 
distributions show distinct peaks at 1.8 ev which are 
insensitive to incident energy, energy density, and 
molecular orientation. Simulations indicate that there may 
be fine structure in the 0-3 ev/sputtered atom range with a 
cascade interaction unique to molecular sputtering. For a 
number of molecular orientations which produce very 
different cascade overlaps, the sputtering yield ratios do 
not correlate with the degree of overlap. The energy 
density of a cascade is not an accurate predictor of the 


sputtering yield ratio. 
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I. INTRODUCTION 


A. OVERVIEW 

Sputtering, the ejection of atoms from a solid by 
energetic ion bombardment, is caused by the collision of 
incident ions with atoms in the surface and near-surface 
layers. 

Sputtering was first observed by Grove [1] in 1852. 
While conducting experiments with electrical discharges in 
a low-pressure gas, he noticed that cathode material was 
deposited on the inner wall of the discharge tube. Thus, 
sputtering gained its reputation as a contaminating 
side-effect of electrical discharge processes. This 
reputation persisted for nearly 100 years and scientific 
interest in the sputtering phenomenon was largely directed 
toward its elimination rather than its understanding. 
However, major scientific achievements in the 20th century, 
such as high voltage, high power radar transmitter tubes, 
focused attention on sputtering. Thereafter, sputtering 
research followed two parallel courses: First, the 
development of a basic theoretical understanding of the 
physics of sputtering, and the acquisition of experimental 
data. Second, exploitation of this knowledge for practical 


purposes. 








Sputtering as a bulk, macroscopic process was 
"understood" and early applied to a limited degree. For 
example, the development of pumps based on ion trapping 
processes and the production of atomically clean surfaces. 
However, very little true understanding of the physics of 
sputtering was achieved through early experimentation or 


theoretical models. 


B. EARLY EXPERIMENTATION 

Experimental results proliferated but, in general, the 
data were plagued by non-reproducibility caused by 
inconsistent experimental conditions [2,3].  Penning and 
Moubis [4] first demonstrated conclusively the effect of 
ambient pressure on sputtering experiments. They also 
coupled a magnetic field to their discharge apparatus to 
reduce back-diffusion of sputtered material to the surface 
[4]. Their experiments produced the first consistently 
reproducible results for energies greater than 500 ev. 

Surface cleanliness, the elimination of oxide layers or 
adsorbed gases, is a critical factor in sputtering 
experiments.  Arifov, et. al., [5] demonstrated in 
controlled experiments the dependence of the sputtering 
yield on surface contamination.  Yonts and Harrison [6] 
presented evidence indicating that surface recontamination 


from background gases is a significant factor in 








quantitative sputtering yield measurements and developed 


criteria for “clean” surface experiments. 


C. DEVELOPMENT OF ANALYTIC THEORY 

Stark [7,8] first extensively discussed the concept of 
an individual sputtering event on an atomic scale and 
developed two competing models: He originated the hot-spot 
model which proposed the ejection of surface material from 
a microscopically small region due to evaporation at high 
local temperatures generated by ion impact [7]. His 
second, the collision model, described the sputtering event 
as a series of binary collisions initiated by a single 
incident 10n. He applied the conservation laws for elastic 
collisions to this model and qualitatively described the 
energy dependence of the yield curve [8]. Stark obtained 
limited agreement with experiments at energies less than 1 
kev. 

Kingdon and Langmuir [9] further developed the binary 
collision momentum and energy model to describe the 
observed desorption of a thorium monolayer from a tungsten 
surface by ion bombardment. 

In 1926, von Hippel [10] formally presented the 
hot-spot theory of sputtering which was further developed 
by Townes [11] in 1944.  Townes' theory correlated well 
with some contemporary experimental results, but the 


agreement was subsequently shown to be fortuitous. 





Keywell [12] developed a theory based upon the 
moderation of neutrons, which may be taken as a precursor 
of the collision cascade concept to follow. Harrison [13] 
developed the first direct statistical theory of sputtering 
which bridged the gap between the Kingdon and Langmuir 
momentum transfer theory and Townes' hot-spot theory. A 
significant feature of the Keywell and Harrison theories is 
the application of probability concepts as expressed by 
collision cross-sections. 

At this point in the evolution of sputtering research, 
the vast majority of experimental data were based upon 
sputtering from polycrystalline surfaces [3]. Sputtering 
was assumed to be isotropic with the angular distribution 
of sputtered material following a cosine distribution. 
Effects due to the underlying crystal structure were not 
considered. 

A major turning point in sputtering theory was the 
experimental demonstration of lattice effects in 
monocrystal sputtering. Wehner [14] showed that sputtering 
from single crystals proceeded in preferred directions 
which correlated with the directions of closest packing in 
the particular crystal. This sputtering anisotropy 
completely contradicted the hot~spot model and pointed to 
momentum transfer as the key ejection mechanism.  Silsbee's 
focusing collision model [15] quickly followed as a partial 


explanation for preferential ejection. In addition, 
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experimental results showed the sputtering yield to be 
Semongly sensitive to the crystallographic direction of the 
incident ion beam. The subsequent discovery of channeling 
[16-20], and the elucidation of assisted focusing and 
replacement collision (previously known as crowdion) models 
[21] showed that even a single momentum transfer theory 


must be an aggregate of many different mechanisms. 


P COLLISION CASCADE THEORY 

Sigmund [22,23] developed an extension of the 
statistical theory which displaced the focusing models. It 
1s based upon a linearized form of the general Boltzmann 
transport equation. Sigmund's model treats an amorphous 
solid with a surface binding energy potential barrier. An 
incident ion strikes a surface atom and creates a collision 
cascade via binary collsions. An atom in the cascade 
sputters if it is headed outward from the surface and 
overcomes the surface potential barrier. In doing so, the 
atom's velocity vector is refracted away from the surface 
normal by the attractive surface binding energy. Signifi- 
Cant agreement with high energy experimental results, that 
is, incident ion energy much greater than the surface 
binding energy, has been achieved. However, it has since 
been shown that the major result of the Sigmund theory fol- 


lows immediately from a dimensional analysis of the system! 


IL 





Thompson [24] and Robinson [25] have taken similar 
approaches while accounting for focusons and channeling, 
respectively, and have obtained comparable results. 

The primary utility of the Sigmund and Thompson 
theories has been to provide experimentalists with simple 
formulae to apply to a wide variety of sputtering 
experiments. However, these theories fail to adequately 
predict sputtering yields at low energies [22,23,26]. The 
high-energy, binary collision model for the collision 
cascade does not account for the predominance of low-energy 
sputtered atoms from near-surface layers. 

In general, the collision cascade can only be 
reasonably understood as a many-body, or multiple 
interaction (MI), process. However, the MI model cannot be 


solved analytically. 


Pe COMPUTER SIMULATION MODELS 

Gibson, Goland, Milgram and Vineyard (GGMV) [27] first 
demonstrated the power and feasibility of the use of 
computer simulation techniques to model radiation damage 
experiments. Their computer program builds a finite, 
three-dimensional monocrystal lattice of copper, 
accelerates a lattice atom in a specific direction with 
specific energy, and solves the resultant MI cascade by 


applying numerical methods to Newton's equations of motion. 
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In this particular case, GGMV studied the generation of 
lattice vacancies and interstitials. 

Robinson and Oen [16, 17] applied GGMV's prototype to 
ion range studies and demonstrated ion channeling, which 
was subsequently verified experimentally [18, 20]. 

Harrison expanded upon GGMV's concept and designed a 
Simulation model to specifically study sputtering [28-32] 
with the flexibility to simulate a wide variety of 
experimental conditions. Harrison's early results showed 
that neither Focusons nor Sigmund's or Thompson's models 
accounted for the majority of sputtered atoms (see also the 
work of Lehmann and Sigmund [33]). Simulations show that 
the vast majority of sputtered atoms originate in the 
surface layer, and the results agree reasonably well with 
experimental yield data. Simulations augment laboratory 
results and allow the isolation and manipulation of 


variables not readily accessible in the lab [34, 35]. 


Ur EIITERATOMIC POTENTIAL FUNCTIONS 

A primary thread of uncertainty that runs through each 
sputtering model and theory is the form of the interatomic 
Berential function £3, 35, 361. No single interatomic 
potential has yet been deduced which fits every sputtering 
Sputtering experiment. Simulations are particularly well- 


Suited for isolating, modifying, and observing the 


13 





sensitivity of the sputtering event to the nature of the 


Nnterstomie potential (35, 37, 38]. 


G. NON-LINEAR CASCADE EFFECTS 

Recent experimental results indicated that the analytic 
sputtering theories do not adequately predict the 
sputtering yield for heavy ions bombarding heavy targets 
[39-41]. A spike theory has emerged which states that 
linear cascade theory breaks down (becomes non-linear) when 
the average energy per atom in the cascade approaches the 
surface binding energy of the target material [42]. The 
spike model emphasizes random ejection effects from the 
region of high energy density. Kelly [43] proposes a 
thermal sputtering theory, a Renaissance of the hot-spot 
theory, to explain the non-linear effects in the context of 
the spike. 

These theories reinforce the multiplicity of sputtering 
phenomena, but continue to fall short of providing a 


comprehensive analytic model. 


H. MOLECULAR SPUTTERING 

Bombardment by molecular ions is well-suited for 
investigating non-linear effects since the yield from a 
Single ion can be compared to the yield from a molecular 
ion of the same element with the same energy per ion. Few 


published results of molecular sputtering exist [39-41, 44]. 
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The available results do demonstrate non-linear effects but 
do not confirm the abnormally high yield enhancement 
addressed by the spike model [42]. 

It is not yet clear whether the individual ion cascades 
resulting from the molecule's dissociation upon impact 
overlap to significantly increase cascade energy density, 
or whether molecular orientation affects the sputtering 


mechanisms. 
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PETA OBJECIIVES 


Few results exist for molecular sputtering. Bader, et. 
al, [44] conducted low-energy (< 8 kev) experiments with 
N* and NÀ on several metals. Their primary 
objectives were to study the effects of the upper 
atmosphere on satellites and other space vehicles and to 
investigate erosion rates within ion propulsion drive 
systems. Recent experiments in non-linear cascade effects 
use high-energy (100-200 kev/ion) polyatomic molecules as a 
convenient method of increasing primary energy density in 
near-surface layers of polycrystalline targets [39-42]. No 
other experiments have been conducted to study the effects 
of molecular sputtering, per se. 

The primary objective of this thesis is to study the 
sputtering of a copper monocrystal by low-energy (0.5 - 5.0 
kev/ion) diatomic molecules using a computer simulation 
model. 

More specifically, the following areas will be 
studied: 

l. Comparison of energy and angular distributions of 
sputtered atoms, and absolute yields for molecular 
(03) and ion (Ot) bombardment. 

2! Effects of molecular orientation upon the 
spüttering yield. 
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3. Comparison of 04/20+ yield ratios with 
published results for N3/2N* [44]. 

fC ecomoacison of 05/20" and Ar s/2Ar” 
yield ratios for indication of yield enhancement in 
accordance with non-linear cascade theory [39-42]. 

5. Effects of mass ratio (M,/M>) and ion size 
(ion-atom potential) upon yield ratios. 

6. Differences in incident ion scattering for 
molecular and ion bombardment. 


7. Cascade overlap. 


L7 





IDE THE TMODEL 


A. BACKGROUND 

QRAD, a full-lattice simulation computer program, is 
designed to model the dissipation of an incident ion's 
momentum, using classical mechanics, in a solid monocrystal 
[28-32]. See references [34] and [35] for detailed 
discussions of the theory and principles of simulation 
physics as applied to sputtering. This thesis is based 
upon the program QMOL, a slightly modified version of QRAD, 
which allows bombardment by diatomic, homonuclear 
molecules. The remainder of this section will be presented 
in terms of QRAD. The unique differences of QMOL will be 
discussed in following sections. 

A note on termionlogy is instructive. Ion refers to a 
Single atom incident upon the target. Molecule refers toa 
diatomic, homonuclear species of the same element as the 
ion incident upon the target. The molecule always 
dissociates upon impact. Immediately following impact, 
each fragment of the molecule is an ion. Atom refers to an 
atom of the monocrystal. A sput is a sputtered atom. 

QRAD simulates a wide variety of experimental 
conditions. Each of the common metallic unit cells can be 


built and duplicated to form a monocrystal of up to 2500 
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atoms with a specific crystal plane exposed to the incident 
ion. Additionally, QRAD can model alloys, stepped 
surfaces, and adsorbed surface layers. 

The ion, of given mass, is shot at the surface by 
assigning to it velocity components consistent with its 
angle of incidence and energy, and positioning it 
immediately above a specific impact point. Hamilton's 
equations of motion are solved numerically by an average 
force method [45] for the resultant MI collision cascade, 
or trajectory. 

The trajectory advances by timesteps. During each 
timestep, the positions and velocities of the ion and atoms 
develop simultaneously. Each ion and atom is "looked at" 
and its motion and forces analyzed at the end of each time 
step. The "look" each receives depends upon its total 
energy and the total energy of its nearest neighbors. A 
hierarchy of three mini-logics corresponding to the depth 
of the "look" significantly reduces computer time compared 
With earlier generations of QRAD. The program flags all 
ions and atoms which rise above the surface. The 
trajectory terminates when the energy of the most energetic 
atom decreases to a specified value below which the 
probability of additional atoms rising above the surface is 
very small. 

Each flagged ion or atom is tested by the program to 


determine whether it has sufficient outward momentum to 
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completely escape the surface into the "vacuum". Atoms and 
ions which pass this test are listed as sputtered and 
scattered, respectively. The program records the final 
positions and velocities of the sputs and scattered ions 
and restores the system so that another trajectory can be 
executed. Individual trajectories may be studied in detail 
to follow the sputtering process or identify the occurrence 
of other phenomena such as channeling. 

Global results which may be compared to experimental 
data are obtained by averaging the results over N 
trajectories, where N is typically 80-100. The sample of 
trajectories is based upon N impact points uniformly 
distributed within an impact area chosen such that it is a 
representative symmetry area of the surface. The N 
trajectories determine sput yields, energy and angular 
Se pDutrons and the spatial distribution of yield over 
the impact area. Multimer formation is determined from 
additional sput analysis. Dimers and trimers are 
identified and listed with their energy and angular 


distributions. 


Be INPUT PARAMETERS 
Molecular data input parameters are listed and defined 
below so that their meanings will be clear in following 


sections. 
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l. Bullet parameter card 


usum ion (2) 


LSS = 0/3 


THEM 


PHIM 


SEP 


IRANF = 0/1 


ISEED 


2 Reconstruction card 
IRECON = 0/1 


ITRI 


ZE 


chemical symbols of 


constituent atoms 


tlag: 
O = single incident ion 
3 = incident molecule 


8m, polar angle of 
ion (2) with respect to 


ion (D) 


$m, azimuthal angle of 
ion (2) with respect to 


Seat 


p, interatomic separation 


in angstroms. 
flag for random molecule 


seed value for random 


no. generator 
flåg for reconstruction 


of trajectory 


no. of trajectory to be 


reconstructed 








p. TMEEROCRYSTAL 

ARSG0O=atom copper (FCC) monocrystalewismused for each 
Simulation. The crystal is "cut" to expose the (111) plane 
to the incident molecule (Fig. 1). The crystal contains 4 
layers with 90 atoms per layer and is designated (15x4x12): 
15 planes wide, 4 planes thick, and 12 planes long.  RYBL 
is the minimum positive distance an atom must rise above a 
surface before it is considered ejected. Note that the 
first numbered atom in the crystal is labelled atom 3 (Fig. 
2). Atoms l and 2 comprise the molecule. When comparing 
OMOL and ORAD trajectories, ATOM(x)omor = ATOM(x-1l)opap- 

Containment, the ability of the target lattice to 
contain the entire collision cascade, is discussed 
elsewhere in detail [34]. The size of this microcrystal is 
a trade-off between asymptotic containment and reasonable 
computer run-time. A series of tests have determined that 
absolute yield containment requires approximately 150 atoms 
per layer at 1.0 kev and about 250 atoms per layer at 5.0 
kev. Absolute containment of an event is paramount for an 
analytic description of the sputtering yield. However, it 
is less significant when studying mechanisms or observables 
such as yield ratios, and energy and angular distributions. 
Experience shows [28,-32, 34] that this lattice provides 
adequate containment for these observables at 0.50 to 5.0 


kev per ion. 
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Pa ORLENTATION OF THE MOLECULE 

Figure 3 shows the orientation of the molecule with 
respect to the crystal axes. Note that the -y axis is the 
outward normal for the (111) plane. 

The molecule has 3 degrees of freedom. Angles 9; and 
$1 define the incident direction of the molecule. 0m and 


gm, the polar and azimuthal angles respectively, and the 


interatomic separation,p , define the position of ion(2) 
relative to ion(1). Ion(2) is assigned the same velocity 
components as ion(l). Vibrational and rotational degrees 


of freedom are not included. 

The ranges of the angles are 0 < 9... 9; < 1/2 
and O < $m, di < 27. Therefore, ion(2) is always trailing 
ion (1) prior to impact. The molecule can assume any 
orientation with respect to the surface because ion(1) and 
ion(2) are identical. This procedure would not be suitable 
for an AB-type molecule such as CO. Each simulation was 
executed with the molecule normally incident upon the 
surface. 0m and ọm were set at specific values for some 
runs, but at random values for others to be more similar to 
experimental conditions. Some effects of beam polarization 
can be inferred from trajectories in which 0m and dm are 
set at specific values : (1) transversely polarized 
molecule, 9m = 7/2, O < dm < 27, (2) longitudinally 


polarized molecule, ðm = m = O. 
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Random values of @m and dm are generated as follows: 
When IRANF is set, the program calls the library random 
number generator LRND which uses the value of ISEED to 
generate 5 random numbers. Four are used to calculate oii 
and $54, by von Neumann's method: 


A aR R 
n 


en = gi a? 
2 2 
(RI + R5) 

QE. Supe 2 R3 R4 
2 
(R3 + R4) 


where O < Sm, Om < 1/2. The fifth sets 0 < dm « "T or 
™< Om < 27. LRND is called for each subsequent trajectory 
and uses its internally updated value of ISEED to generate 
the next set of random numbers. 

The initial position of ion(2) with respect to the 


impact point is calculated as follows: 


DPX = pp sin 8m cos om 
DPY = -p cosd m 
PEZ- Sino m Sino m. 


The exact impact point of ion(2) is unknown prior to impact 
unless the molecule is transversely polarized. In all 
cases, ion(2)'s motion is perturbed by interaction with 
ion(1) at impact.  Ion(2)'s impact point can be determined 
by examining the first 5-10 timesteps in the printed output 
of a trajectory reconstruction. 

A random molecule trajectory is not reproducible unless 


the exact values of 9m, óm, and ISEED are repeated from the 
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Original trajectory. A specific trajectory is 
reconstructed as follows: When IRECON is set, the program 
calls LRND and executes a LRND loop ITRJ-times, where ITRJ 
1s the number of the trajectory to be reconstructed. The 
program now holds the same set of 5 random numbers for 

aN CoV RR as an the original run and calculates 9. 


and $4 as before. 


E; IMPACT AREA 

The impact area for each simulation is a rectangle 
containing a grid of 104 uniformly distributed impact 
points (Fig. 4) and is a representative symmetry area of 
the entire (111) surface (Fig. 2). The impact area was 
designed for bombarding a stepped Ni (111) surface [46] and 
was adapted to this research to provide an extensive 
sampling area. 

Note that the reference target point is centered on 
atom 36 (for QMOL) and is designated (RBX, RBZ) = (6, 4) in 
units of planes as measured from the origin at the center 
of atom 3. Likewise, the impact points, (SPX, SPZ), are in 
units of planes and are labelled in figure 4 with their 
corresponding trajectory numbers. Figure 5, drawn to 
scale, shows the complete grid of impact points 


superimposed on the surface with atom 44 centered in the 
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ELesNote tae positions Of the atoms in the second 
layer. 

Independent data sets can be generated from the 
104-point data set by simple displacements from the 
original RBX/RBZ. Figures 6 and 7 show two sets of RBX/RBZ 
displacements which can generate 1-12 displaced data sets. 
The utility of the displaced data sets is to significantly 
increase the sput population for a particular simulation 
without significantly changing the position of the impact 
area. The increased populations reduce statistical 


flucuations in energy and angular distributions. 


Pee NE RGLES 

Recall that Bader, et. al., [44] conducted low-energy 
(< 8 kev/ion) molecular sputtering experiments, and that 
investigations of non-linear cascade effects have utilized 
high-energy (100-200 kev/ion) molecules. The simulations 
presented in this thesis were conducted at 0.50, 1.25, 

EE anes .00 key /ion. 

This is the first investigation of molecular sputtering 
using the latest-generation simulation model. An extensive 
data base of single-ion simulation sputtering experiments 
exists in this energy range. A molecular sputtering data 
base of similar scope will provide additional insight into 


the analytic discrepancies of low-energy sputtering. 


28 








The results of these simulations can be compared with 
mae Experimental results of Bader, et. al., [44]. 

Andersen and Bay suggest that non-linear yield 
enhancement is a slowly decreasing function of energy [39]. 
The trends of the 04/207 yield ratios at these low 
energies will be instructive. 

Note that the heat of sublimation of copper is 3.53 
9eV/atom and the surface binding energy is 2.45 
Sato for the potential functions set into the model 
[37]. Assume that the incident ions (QRAD) and molecules 
(QMOL) deposit all their energy into the crystal. The 
resulting energy per atom values will bracket the heat of 
sublimation and surface binding energy values. Systematic 
effects of the relative magnitudes of the energy density, 
and the surface binding energy and heat of sublimation in 


the context of the spike model [42] can be studied. 


GT POTENTIAL FUNCTIONS 
1. Atom-Atom potential functions 
The program uses the Pot-II potential function, a 
compound function which provides a repulsive Born-Mayer 
(Gibson-2) "wall joined smoothly to an attractive Morse 
"well" by a cubic spline at the points Ra = 0.83 LU and 


[p XI. 9NIU. See figure 8. 
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The repulsive wall dominates the collision cascade 
Gymamics. Attractive forces have little effect upon the 
dynamics but dominate the atom's ability to escape the 
surface. 

The potential function well depth reproduces the 
solid's heat of sublimation (total energy) by summing ever 
all possible pairwise interactions. The Cu-Cu potential 
includes only nearest neighbor and next-nearest neighbor 
Beacons because JE 1S truncated at a radius R, = 
2.40 LU. The well also reproduces the solid's static 
surface binding energy (SBE) at its pre-trajectory level. 
The SBE is the attraction a sput must overcome. The values 
of the heat of sublimation and SBE are 3.53 ev/atom and 
2.45 ev/atom, respectively. 

See reference [38] for a detailed discussion of the 
atom-atom potential function. 

Ze Ton Atom Potential Function 

The Cu/O* ion-atom potential function is a 
compound functon consisting of a 0.7-Moliere wall splined 
to a Morse well at Ra = 0.44 LU and Rp = 0955 50. 

Refer to reference [37] for a detailed discussion 
of the ion-atom potential function. 

eon Llon Potential Function 
The oxygen molecule potential function is a Morse 


potential truncated at R, = 2.40 LU (Fig. 8). 
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IV. DISCUSSION AND RESULTS 


A. ENERGUUDISTRIBUTTONS 

The characteristics of the sputtered atom energy 
wee butions for the incident energies of 0.5, 1.25, 2.5, 
and 5.0 kev/ion are compared for the following simulations: 
(1) QRAD, (2) random molecule, QMOL(R), and (3) longitu- 
dinally polarized molecule, QMOL(LP). These distributions 
are equivalent to those which would be obtained from an 
ideal experiment which bombards a clean, monocrystalline 
(111) copper surface in a perfect vacuum, and collects and 
analyzes each sputtered atom. 

1. Normalized Energy Distributions, 0-20 ev/atom 

Figures 9-12 are normalized energy distribution 

curves comparing QRAD, QMOL(R), and QMOL(LP) for each 
incident energy and a sputtered atom energy range of 0-20 
ev/atom. These curves were obtained from total sputtered 
atom populations of 700-1800 atoms and are plotted as raw 
data. Individual sputtered atom populations are indicated 
with each distribution. These small populations produce 
distributions with many structural details which are not 
statistically significant. The smooth curves fitted to the 
data have no physical significance for any of the 
distributions taken alone. However, when viewed in a 


series, certain structural trends may be elucidated for 
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specific simulations at different incident energies or 
different simulations at a specific energy. Several 
features of the distributions are strikingly apparent. The 
maxima occur at precisely the same energy, approximately 
1.8 ev/atom. The only exception is the maximum for QRAD at 
0.5 kev/ion. However, its maximum at approximately 2.2 
ev/atom is within the statistical uncertainty of the 
simulation and the difference represents a flucuation of 
8-9 sputtered atoms. 

A shoulder appears in each distribution between 2 

and 5 ev/atom. The shoulders for QRAD and QMOL(LP) 
decrease in energy and increase in fraction of total yield 
as incident energy increases. The shoulders for QMOL(R) 
remain nearly constant as incident energy increases: The 
widths of the shoulders are 1-1.5 ev/atom and the yield 
fractions are 0.040-0.054. This is equivalent to a yield 
flucuation of 41-95 sputtered atoms. Secondary shoulders 
may be present in the 6-10 ev/atom range as well as valid 
structure at higher energies, but the populations are 
insufficient to define additional features. 
2. Normalized Energy Distributions, 0-3 ev /atom 

It is instructive to study the normalized energy 
dusctsrsubrons for the sputtered atom energy range 
containing the peaks seen in the preceding section. 
Figures 13-16 are distributions for the energy range 0-3 


ev/atom. These distributions are plotted from a "sliding-4" 
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average of the raw data as follows: raw data points 1-4 
are summed, followed by points 2-5, 3-6, 4-7, etc. The set 
ef averaged points is normalized and plotted as the 
fraction of total yield vs. energy per sputtered atom. As 
before, the smooth curves fitted to the data have no 
EudrIvrdual physical significance. 

The relative locations of the apparent maxima of 
the distributions are consistently skewed toward lower 
energies as indicated by lines A in the figures. 
recante shifts of the distributions for a particular 
Simulaton as a function of incident energy are not seen. 
In fact, the constant energy peaks of figures 9-12 suggest 
Bueno significant shifts should occur in the 0-3 ev/atom 
range. 

Several apparent peaks are seen in the distribu- 
tions. Consider the QRAD distributions. The peaks are 
fitted credibly because of their consistency as the 
incident ion energy increases. The peak vanishes at 5.0 
kev/ion, but note that the total sputtered atom population 
at 5.0 kev/ion is only half of the populations at the lower 
incident energies. The 5.0 kev/ion QRAD simulation 
consists of 1352 trajectories. Nearly 3000 trajectories 
are necessary to produce a sputtered atom population 
sufficiently large to verify a peak. No comparable peak 


behavior is seen for QMOL(R) or QMOL(LP). 
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Several low-energy shoulders may be real and 
verifiable for significantly increased populations. 

Assume that these structural features are real. 
Individual sets of trajectories isolating the atoms 
comprising the features can be studied to: (1) identify 
the sputtering mechanism generating the QRAD peaks and 
shoulders, and (2) elucidate the molecular effects 
responsible for breaking up or reinforcing the mechanisms. 

3% Incident Energy Dependence of the Normalized Energy 
Distributions 

Figures 17-22 show the ion/molecule energy 
dependence of the normalized energy distribution functions 
for each simulation. The distributions are identical to 
those previously studied and are displaced vertically along 
an incident energy scale. Fraction of total yield cannot 
be read from this scale but an increment of fractional 
yield is indicated for each figure. 

The 0-20 ev/atom distributions, figures 17-19, 
clearly show the alignment of the peaks and the shifts in 
shoulder positions. The QRAD shoulders shift linearly to 
lower energy and higher fraction as incident energy 
increases. The QMOL(R) shoulders remain constant in energy 
with Wen icant shite in y te ld thaetton. e O MOL (EE) 
shoulders show a non-linear shift to lower energy and 


higher fraction. 
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Pace cuz W=22 sare thewo-3 ey/atom distributions. 
The QRAD peaks remain in the vicinity of 2 ev/atom. The 
prospective location of a peak in the 5.0 kev/ion 
distribution, discussed in the preceding section, is clear 
in figure 20. No significant trends in fine structure are 
evident in the QMOL(R) or QMOL(LP) distributions. Real 
shifts in the distributions for either simulation, as a 
function of incident energy, are not supported. 

a lative Energy Distribution Functions 

Figures 23-25 are the cumulative energy distribu- 
HEcpDESunetlons (CDF) which plot the fraction of sputtered 
atoms with energy greater than the indicated energy as a 
function of sputtered atom energy for the range 0-200 
ev/atom. The figures compare the simulations as follows: 
eee 2s, ORAD vs. QMOL(R), (2) Fig. 24, QRAD vs. 
eromımE), and (3) Fig. 25, QMOL(R) vs. OMOL(LP). The 
incident ion energy is indicated with each curve. Note 
that for a specific incident ion energy the lower of the 
two curves shows a shift to the low energy region of the 
range since each plotted point represents the fraction 
greater than the energy indicated. 

Figure 25 shows the only consistent trend for the 
distributions. The QMOL(R) sputtered atom populations are 


shifted to lower energies than those for QMOL(LP). 
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Figure 23 will be particularly instructive in 
eee-elating the angular distribution effects discussed in 


the following section. 


B. ANGULAR DISTRIBUTIONS 

Figures 26-29 are the angular distributions of 
sputtered atoms (spot patterns) for QRAD and QMOL(R). 
These figures are equivalent to spot patterns obtained 
experimentally by collecting sputtered deposits on flat 
glass plates. 

The most obvious feature of the spot patterns is the 
six-fold symmetry which has been shown experimentally [14] 
to be a characteristic of preferential ejection from the 
(111) surface. The diffuse distribution between the 
six-fold axes is caused by low-energy (<20 ev) sputtered 
atoms ejected late in the trajectories when extensive 
surface disruption has broken down the preferred ejection 
mechanisms [47]. 

The density of the spot patterns is directly related to 
the sputtered atom population. However, there are no 
consistent trends for QRAD vs. OMOL(R) as a function of 
incident energy. The QMOL(R) spot patterns are not always 
"darker" or more diffuse than the corresponding QRAD spot 
patterns. 

The relative densities of the spot patterns for QRAD 


and QMOL(R) for each incident energy can be correlated to 
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ue the CDF for QRAD vs. QMOL(R). Consider the 
0-20 ev/atom energy range for each incident energy. 
IO kev/ 1on 
me ix fold axes for ORAD in figure 26a are very 
well defined but the spot pattern for QMOL(R), figure 26b, 
is heavily filled between the axes. Figure 23a shows that 
the CDF for QMOL(R) is shifted more toward lower energies 
than for QRAD. The higher fraction of low-energy sputtered 
atoms is responsible for the heavy, diffuse distribution 
between the axes. 
pee. 2 > key /ion 
The CDF for 1.25 kev/ion incident energy, figure 
Ps exactly opposite that for 0.5 kev/ion but the 
differences in fraction are not as great. This is 
reflected in the spot patterns, figures 27a and 27b. Each 
is very dark but the QRAD pattern now shows slightly more 
filling between the axes than the QMOL(R) pattern. 
Eee) «kev /ion 
The CDF and spot patterns for 2.50 kev/ion incident 
energy, figures 23c and 28, respectively, follow the same 
trend as those for 1.25 kev/ion incident energy. The 
fractional differences in the CDR are smaller but the 
higher QRAD sputtered atom population produces a darker 
spot pattern ( the percentage difference in sputtered atom 
population between QRAD and QMOL(R) is greatest at this 


incident energy). 
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SS Oe 9n 
The CDF curves for 5.0 kev/ion incident energy, 
Zeure 23d, are virtually identical, and the spot patterns 
for this energy, figures 29a and 29b, are nearly 


indistinguishable. 


©. SPUTTERING YIELD RATIOS 
Figure 30 compares the experimental yield ratios, for 
the sputtering of polycrystalline copper by N* and 
N*5 [44], and the simulation sputtering yield ratios. 
The sputtering yield ratio is defined as follows: 


ratio = yield from X9T 


Ne o) 


where X is the bombarding element. Two sets of simulations, 
"real" and "artificial", are considered. 
l. Experimental 


The experimental sputtering yield ratios for N* 
and N5* incident upon polycrystalline copper [44] are 
plotted in figure 30a and 30b for the incident energy range 
0-8 kev/ion. A very good straight-line fit is obtained 
which shows a slowly increasing yield ratio with increasing 
incident energy. 

2. "Real" Simulations 
The sputtering yield ratios for the "real" 


simulations, QMOL(R) and QMOL(LP), are plotted in figure 
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30a. The connecting lines have no physical significance. 
Maegeneral, the simulation yield ratios increase much more 
slowly with incident energy than the experimental values. 
The ratios for OMOL(LP) are consistently lower than those 
"or OMOL(R) and both dip considerably at 1.25 kev/ion 
imeident energy. 
EN artificial” Simulations 

Hicmsputtering yield ratios for two "artificial" 
Simulations are plotted in figure 30b. In the first case, 
Ar* and Ar} are substituted for 0* and 0 っ + 
but the Cu-O0 ion-atom potential is retained. Therefore, 
the mass ratios are changed for QRAD and QMOL(R) from 0.251 
ee 0, for 07/05" to 0.628 and 1.257 for 
Ar*/Ar5*, respectively. In the second case, the 
B-potential [37] is substituted for the Cu-0 potential and 
the oxygen mass is retained. The ratios for Ar3/Ar* 
are lower than the corresponding "real" QMOL(R) ratios in 
figure 30a and an identical dip occurs at 1.25 kev/ion is 


almost totally absent. 


Pees beers OF MOLECULAR ORIENTATION 

The effects of molecular orientation upon the 
Sputtering yield ratio are examined for the following: (1) 
longitudinally polarized molecule, QMOL(LP), (2) 
transversely polarized molecule (8m = 90 deg.), QMOL(TP) 


and (3) polar angle 9m = 45 deg., QMOL(45). Figures 3la 
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and 3lb show the Sputtcaingmyreld ratio as a function of 
azimuthal angle, fm» at 1.25 kev/ion incident energy for 
QMOL(TP) and QMOL(45). Lines of constant yield ratio are 
indicated for QMOL(LP) at 1.25 kev/ion incident energy and 
the average QMOL(LP) yield ratio for the four incident 
energies. Note that azimuthal angle is not defined for a 
longitudinally polarized molecule. 

1. General 

On the average for a random molecular orientation, 
the two dissociated ions will not collide with the surface 
simultaneously. When the molecule is transversely polar- 
ized, the time separation between collisions is minimized 
and the separation between impact points is the greatest. 
The degree of overlap of the individual collision cascades 
will be determined by the separation in the diatomic 
molecule, 1.22 A for O2. When the molecule is 
longitudinally polarized, the time separation is maximized 
and the impact points are very nearly identical. 

The longitudinally polarized molecule produces the 
highest energy density in the collision cascade since the 
energy from both ions is deposited in the same region of 
the lattice. If energy density were the controlling 
factor, the longitudinally polarized molecule should 
produce the highest yield ratios. 

Mentra Mng ion in the longitudinally polarized 


molecule is colliding with a surface already disrupted by 
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the leading ion and, in general, will not "see" a surface 
ES tonto tone Impact. In adition, the trailing ion 
overtakes the leading ion and is deflected slightly by the 
leading ion's collision with the surface. These factors 
may result in increased penetration depths and a higher 
incidence of channeling for bombardment by longitudinally 
polarized molecules. 
2. Yield Ratios 
Mie yteola ratios for OMOL(LP) and O@MOL(45) at 1.25 
kev/ion incident energy are plotted in figure 31, with 
OMOL(LP) yield ratios indicated, as a function of the 
azimuthal angle S The ratios for both cases show a good 
linear fit and are nearly identical. The QMOL(LP) values 
are not significantly different from either QMOL(TP) or 
QMOL (45). 
3. Scattered ions 
The following table summarizes and compares total 
yields and the number of scattered ions for QRAD, QMOL(R) 
and QMOL(LP): 


no. 


traj total yield scattered ions 




























Om) kev/ion QRAD 936 1665 104 
QMOL(R) 520 LITO 223 
QMOL (LP) 520 1657 138 
MaS Lon QRAD 936 1653 3g 
QMOL(R) 520 1660 67 
QMOL (LP) 520 1437 42 
ISO WeYe QRAD 1352 1439 20 
QMOL(R) 520 1181 29 

QMOL (LP) 520 1069 
8 


5.0 kev/ion QRAD 352 927 
QMOL (R) 520 
QMOL (LP) 520 
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eee OLHER MOLECULAR EFFECTS 
l. Yield per Layer 

Detailed analyses of single-ion sputtering 
trajectories have shown that the majority of sputtered 
atoms come from sputtering mechanisms which are completely 
confined to the surface and near-surface layers [29]. 
Consequently, the majority of sputtered atoms come from the 
surface layer. 

The sputtering yields for each layer of the 
microcrystal are listed below, for each simulation and 
incident energy, with the fractional yields from layer 2 
(Lo2/N) and all subsurface layers ((L> * L4 * L4)/N) 


ined cated < 


Ges Key /ion N ee ie oy eo N N 


1635/28/2/0 
1708/53/7/2 





0.016 0.018 
029 0035 


QRAD 1665 
QMOL(R) 1779 





QMOL(LP) 1657 1989/547/3 O 0317 or od 
1653 1626/23/2/2 -013 0.016 
OMOL ( R ) 1660 1614/31/15/0| 0.018 02.027 
QMOL(LP) 1437 1397/35/5/0 .024 07007 
1439 1420/16/3/0 Sog 3 
OMOL ( R ) 1181 12197, 317370 .026 0.028 
QMOL(LP) 1069 TOS” SSI O 2032 0.033 
3.0 kev/ion 
QRAD 927 909/17/1/0 0.018 0.019 
OMOL(R) 760 oA dA 0.027 0.030 
OMOL(LP) 746 718/26/2/0 0.034 05037 
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Esteorhat For molecular sputtering, the vast majority of 
Senec ered atoms still originate in the surface layer. 
However, the relative yields from subsurface layers are 
much greater for molecular than for single-ion sputtering. 
2. Normalized Sputtering Frequency 

Figures 32-39 are the surface layer sputtering 
frequency-location diagrams for QRAD and QMOL(R) at each 
incident energy. Each diagram is normalized to the most 
frequently sputtered atom for the simulation: a circle 
containing ten rings indicates the most frequently 
sputtered atoms(s). A circle containing zero rings is an 
atom that never sputtered and one ring indicates an atom 
atom that sputtered at least once in the entire set of 
trajectories. 

Atom 43 (44), for QRAD (QMOL), is the primary 


knock-on atom (PKA) in at least 46% of the trajectories 


(see figure 5). It is always one of the least sputtered 
atoms. In fact, its frequency decreases as incident energy 
increases. Note that the PKA sputtering frequencies are 


always equal for QRAD and QMOL(R) at each incident energy, 
with the exception of 2.50 kev/ion: at this energy the 
frequencies differ by less than 0.1. 

At 0.5 kev/ion, the most frequently sputtered atoms 
are both nearest and next-nearest neighbors, relative to 
the PKA. But at 1.25, 2.5, and 5.0 kev/ion the most 


frequently sputtered atoms are almost exclusively next to 
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the most frequently sputtered atoms are triangular or 


diamond-shaped. 


a MODEL 

Computer run-times for QMOL are, on the average, nearly 
twice those for QRAD at the same energy per incident ion. 
This is reasonable since a QMOL simulation is essentially a 
superposition of two single-ion collision cascades.  Run- 
time increases as incident energy decreases because of ion 
size effects. Several QMOL simulations, particularly at 
m Oen iion, had to be split into two separate runs to 
avoid exceeding the computer system's job time constraints. 

The reproducibility of a QMOL trajectory is extremely 
sensitive to the precision of ?m and ?m. The IBM 360 
system is designed to retain 7 significant digits, in 
Single-precision, and this is the precision of the random 
angles 9m and ŷm. An attempt to reconstruct a trajectory 
by manually entering ®m and $m to 3 significant digits 
produced Only 2 of the original 10 sputtered atoms for the 
trajectory. Subsequent entries of 4 and 5 significant 
digits produced yields of 4 and 7 sputtered atoms, 
respectively. The reconstruction loop reproduces a 


trajectory exactly. 
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V. CONCLUSIONS 


B THE MODEL 

This research reinforces the feasibility and 
E bITty of computer simulated sputtering. The single- 
ion sputtering model was easily modified to simulate 
homonuclear, diatomic molecular sputtering within the 
original constraints and basic physical assumptions. No 
additional approximations or rethinking of the physics 


inherent in the basic model were necessary. 


Br SEUN BERED ATOM NORMALIZED ENERGY DISTRIBUTIONS 

Recall the alignment of the peaks in the 0-20 
Br eanttered atom distributions. The peaks are clearly not 
an incident energy, energy density, or molecular orienta- 
tion effect. Preliminary results of other simulations, not 
included here, show that the peaks are insensitive to the 
mass ratio, but do shift to lower energies when the 
ion-atom potential function is changed from the Cu-0* 
compound potential to the B-potential. The shift is nearly 
identical for single-ion and molecular sputtering. 

The 0-3 ev/sputtered atom energy distributions provide 
the greatest insight into the effects of the molecule upon 
the sputtering process. The ORAD peaks for 0.50, 1.25, and 
2.50 kev/ion incident energy may be real. A mechanism 


unique to the QMOL cascades is clearly breaking up the 
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peaks and in one case (QMOL(R) at 2.50 kev/ion) may be 
reinforcing the peak. The source of this molecular 
sputtering mechanism may be linked to the substantially 
increased fraction of sputtered atoms from subsurface 


layers for the QMOL simulations. 


eo) LLB RATIOS 

The simulation sputtering yield ratios agree reasonably 
well with experimental data. Two primary factors may cause 
the differences in yield ratios: (1) the difference in ion- 
atom potential functions, and (2) the polycrystalline 
target used in the laboratory experiments. Mass ratio is 
probably not significant since the mass ratios for N and O 
Pe Fersbyzonly 12%. 

The majority of the yield ratios are less than one and, 
in general, the overall trend for the yield ratios is a 
slow increase with increasing incident energy. Energy 
density is not an accurate predictor of the sputtering 
peld ratio. In fact, the yield ratios for the 
longitudinally polarized molecule, the configuration 
producing the highest energy density in the lattice, are 
all lower than the ratios for the corresponding random 
molecule simulations. 

The anomalous dip in the sputtering yield ratios at 


1.25 kev/ion may be an artifact of the ion-atom potential 
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function since it vanishes in the B-potential "artificial" 
Simulation. It is interesting to note that 1.25 kev 
averaged over the 360 - atom microcrystal is only 0.06 
ev/atom less than the heat of sublimation of the 


microcrystal modelled by the Cu-Cu potential function. 


D. MOLECULAR ORIENTATION 

spúrBtering yield ratio is strikingly insensitive to 
molecular orientation at 1.25 kev/ion incident energy. 
Although QMOL(TP) and QMOL(45) yield ratios are not 
available for the other incident energies, no significant 
difference is expected since the QMOL(R) and QMOL(LP) 
ratios are close at those energies. 

Molecular beam polarizability may have a definite 
application for ion implantation. The QMOL(LP) simulations 
produced substantially less scattered incident ions at 


NEP and 5.0 kev/ion- 
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VI. RECOMMENDATIONS 


A. THE MODEL 

QMOL can easily be modified to simulate sputtering by 
an AB-type molecule as follows. Change the random number 
generator loop to produce six random numbers instead of 
Enc Use the sixth random number to set O < 6m < m or mv < 
0m < 27 in the same way that the fifth random number sets 
the value of $m. The atom that is in the lower hemisphere 
of the molecule's locus of orientation becomes ion(1). The 
program then sets up the molecule at the impact point. In 
principle, the program can be modified to simulate 
ERaucertersng by triatomic or higher-order polyatomic 


molecules. 


Ba SPUTTERED ATOM NORMALIZED ENERGY DISTRIBUTIONS 

It is necessary to identify specific atoms comprising 
the QRAD peaks. Then, individual QMOL trajectories common 
to these atoms can be examined in detail to elucidate how a 
peak is dispersed into a much wider energy band. However, 
sputtered atom populations for QRAD at 5.0 kev/ion and all 
QMOL simulations are too sparse and statistically uncertain 
Lom ussedetailed analysis. The 5.0 kev/ion QRAD 
simulation should be run with a new displaced data set 


denned tOo augment the 1352 trajectories already executed 
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and produce a sputtered atom population sufficiently large 
to confirm or negate the appearance of a peak. Likewise, 
each QMOL sputtered atom population should be augmented by 
runs with the displaced data set no. 2 impact points. The 
combination of much larger sputtered atom populations and 
the sliding-4 averaging technique will enable more precise 
correlation between structural features in the 0-3 


ev/sputtered atom energy distributions. 


ere YIEDbD RATIOS 

A systematic study of the differences between the 
simulation and experimental sputtering yield ratios must 
focus upon variations in the ion-atom potential function 
until laboratory data for single-crystal molecular 
sputtering becomes available. Two primary options are 
available: (1) Model the Cu-N* potential while retaining 
the oxygen mass in the simulation, and (2) increase the 
O9 bond length to reduce cascade overlap. 

Individual trajectories should be compared for QMOL(R) 
EDEN os nd 1.25 kev/ion, and QMOL(LP) at 1.25 and 2.50 
kev/ion to bracket the dip at its largest variation. The 
relative positions of the ions for each type of simulation 
can be traced from timestep to timestep in the individual 
trajectories. This may show peculiar differences in 


cascade overlap causing inconsistent yield ratios. 


47 





D., TON=ATOM POTENTIAL FUNGCETON 

The cubic spline in the Cu-OT potential function must 
be adjusted. A slight knee in the curve occurs at approx- 
imately 100 ev potential energy and 0.5 LU separation when 
ehe function is plotted from the program's potential table. 
This shows that the spline is not precisely matching the 
magnitude and slope of the Moliere and Morse potentials at 
Piemecimm@ent values Of Ra = 0.44 LU and Rp = 0.55 LU. 
For an initial iteration change Ra and Rp to 0.40 and 
0.60 LU, respectively, execute the program for one timestep 
to print a new potential table, and replot the ion-atom 
potential function. A change in the spline will not affect 


Baessiymulation results. 
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Fig. 5. Molecule orientation with respect 
to the (111) surface of the target. 
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Fig. 4. Impact area and impact point coordinates. 
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Fig. 5. Impact area on the (111) surface, drawn to 
scale. The central atom is atom 44 (QMOL). 
The dashed circles are second-layer atoms. 
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SHIFT 


0.14286 





a = 0.02564, b = 0.02381 


POINT RBX RBZ 
1 6.0000 4.0000 
2 0202 21 4.02564 
5 Ga sod 5.97436 
4 5.97619 5.97436 
5 Soe 4.02564 


All units are planes. 


Fie- 6. Displaced data set no. 1. 
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SHIFT 


0.14286 





Os So 020-272 

POINT RBX RBZ 
6 5.96428 3.96154 
7 5.96428 4.00000 
8 5.96428 4.03846 
9 6.00000 4.03846 
10 6.03572 4.03846 
11 6.03572 4.00000 
12 6.03572 3.96154 
13 6.00000 3.96154 


An ies are planes. 


Pee Displaced data seti no: 2. 
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Fig. 9. Normalized energy distribution of sputtered atoms, 0-20 ev, 
0.50 kev/ion. 
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Fig. 10. Normalized energy distribution of sputtered atoms, 
0=20 "ev, 1.25 kev/1on. 
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Fig. 11. Normalized energy distribution of sputtered atoms, 
0-20 ev, 2.50 kev/ion. 
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Fig. 12. Normalized energy distribution of sputtered atoms, 
0-20 ev, 5.0 kev/ion. 
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Fig. 13. Normalized energy distribution of sputtered atoms, 
0-3 ev, 0.50 kev/ion. 
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Fig. 14. Normalized energy distribution of sputtered atoms, 
0-3 ev, 1.25 kev/ion. 
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Fig. 18. QMOL(R) normalized energy distributions of sputtered atoms, 
0-20 ev. 
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Fig. 21. QMOL(R) normalized energy distributions of sputtered atoms, 
0-3 ev. 
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Fig. 22. QMOL(LP) normalized energy distributions of sputtered atoms, 
0-3 ev. 
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Fig. 23. Cumulative energy distributions of sputtered atoms. 
( + QRAD, o QMOL(R)) 
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Fig. 25. Cumulative energy distributions of sputtered atoms. 
( + QMOL(R), o QMOL(LP)) 
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